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Abstract. The majority of extrasolar planets discovered to date have sig- 
nificantly eccentric orbits, some if not all of which may have been produced 
through planetary migration. During this process, any planets interior to such 
an orbit would therefore have been susceptible to resonance capture, and hence 
may exhibit measurable orbital period variations. Here we summarize the re- 
sults of our investigation into the possibility of detecting low-mass planets which 
have been captured into the strong 2:1 resonance. Using analytical expressions 
together with simulated data we showed that it is possible to identify the exis- 
tence of a low-mass companion in the internal 2:1 resonance by estimating the 
time-dependant orbital period for piecewise sections of radial velocity data. This 
works as long as the amplitude of modulation of the orbital period is greater 
than its uncertainty which in practice means that the system should not be too 
close to exact resonance. Here we provide simple expressions for the libration 
period and the change in the observed orbital period, these being valid for arbi- 
trary eccentricities and planet masses. They in turn allow one to constrain the 
mass and eccentricity of a companion planet if the orbital period is sufficiently 
modulated. 



1. Resonance theory 

The following summarizes the results of Kennedy & Mardling (in preparation). 
A resonant system is characterized by the libration of one or more resonance 
angles. For a coplanar system the quadrupole 2:1 resonance angle is given by 

(j) = Aj - 2A D + Wi (1) 

where \ and A D are the inner and outer mean longitudes, and W{ and w a are 
the corresponding longitudes of periastron. The maximum possible variation in 
a = Vijvo, the ratio of inner to outer orbital frequencies, is given by (Mardling 
2008) 

Aa = 3 [f( ei ) g{e ) ((m /m*) + 2 2 / 3 (m i /m„))] 1/2 , (2) 

where and e D are the inner and outer orbital eccentricities, mi and m are the 
corresponding planet masses, 

f(e i ) = 3e l -^-^ + 0(eJ) and g(e ) = 1 - \e\ + ±f + 0( e °). (3) 
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For stable systems inside the 2:1 resonance the libration period is approximately 

P hb = 2P /Aa = a[f(e i )]- 1 / 2 P , (4) 

where a depends on observed parameters when rrij <C m a . The amplitude of 
variation of the outer (observed) orbital period depends on the distance from 
exact resonance, 5a, and is given by 



SPo/Po 



1 + (m /mi)a 2/3 1 5a ja ~ -2 2/3 (mj/m G ) 5a, 



(5) 



where 5a = a — 2. If the period and amplitude of variation of the orbital period 
of a known system can be measured, (HJ) and ([5]) can be used to estimate ei, nii 
and 5a. 



2. Data Analysis 

To illustrate our procedure we used the orbital parameters of HD 216770b for 
which m D sini = 0.65M/, M* = O.9M , e Q = 0.37, P Q = 118.45 days and 
Pub — 25P . Using synthetic data for a system with and without a 10 Earth- 
mass planet in the interior 2:1 resonance, we were able to infer the presence of 
the super-Earth as follows. First we determined the range of values of for 
which stable orbits exist using direct three-body integrations. Sampling at the 
rate of around 16 synthetic radial velocity data points per orbit and five orbits 
per segment, we found single-orbit solutions to each segment using a Levenberg- 
Marquardt algorithm. The resulting successive orbital periods are then fit by a 
sinusoid and the period and amplitude compared with ([3D and (0). Using this 
procedure we were able to recover the true libration period and amplitude to 
within 30%. This is illustrated in Figure [TJ 
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Figure 1. Piecewise estimates (filled circles) of the orbital period for a sim- 
ulated HD 216770-type system (a) with and (b) without a 10M® companion 
in the interior 2:1 resonance. The solid curves are from three-body integra- 
tions while the dashed curves are sine curve fits to the points, constrained 
by knowledge of possible values for the libration period. Simulated radial 
velocity data includes instrument and stellar jitter noise. 
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